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The relationship between solar wind induced signal phase fluctuation and solar wind 
columnar electron density has been the subject of intensive analysis during the last two 
decades. In this article, a sizeable volume of 2.3-GHz signal phase fluctuation and 
columnar electron density measurements separately and concurrently inferred from 
Viking spacecraft signals are compared as a function of solar geometry. These data 
demonstrate that signal phase fluctuation and columnar electron density are proportional 
over a very wide span of solar elongation angle. 

A radially dependent electron density model which provides a good fit to the columnar, 
electron density measurements and, when appropriately scaled, to the signal phase 
fluctuation measurements, is given as N^(r) = 2.21 X 10^ + 1.55 X 10^ r~^'^ , where 

r = heliocentric distance in solar radii and NJr) = electron density in cm~^. This model 
is also in good agreement with K-coronameter observations at 2 solar radii (2r^), with 
pulsar time delay measurements at lOr^, and with spacecraft in situ electron density 
measurements at 1 AU. 


I. Introduction 

During the last two decades, solar wind electron density 
and density fluctuation have been extensively investigated and 
analyzed via a wide variety of direct and indirect techniques. 
Although a number of studies have been performed using in 
situ spacecraft measurements at approximately 1 AU (e.g., 
Refs. 1-3), the very interesting regions close to the Sun were 


and continue to be inaccessible to spacecraft in situ measure- 
ments, and hence can only be probed indirectly via effects on 
observed signals which pass through or originate in these 
near-Sun regions. 

In the early to mid 1960s, solar wind electron density 
fluctuation experiments focused on measurement of the 
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angular broadening of compact natural radio sources and subse- 
quent modeling of this parameter as a function of solar elonga- 
tion angle (Refs. 4-6). In the latter 1960s the emphasis shifted 
to measurement of signal intensity fluctuation in the form of 
the “scintillation index m” of natural radio sources (Refs. 
7-11). Later studies (Refs. 12-16) emphasized the determina- 
tion of the solar elongation angle (or radial) dependence of the 
scintillation index. More recently obtained monochromatic 
source spacecraft signal scintillation measurements are given in 
Ref. 17. Very recently, near-Sun spacecraft monochromatic 
signal spectral broadening has been studied (Refs. 18-20). 

In the early attempts to formulate a theoretical basis for 
natural source intensity scintillation, it was necessary to con- 
sider signal phase fluctuation induced by electron density 
fluctuation. Studies which investigated the geometric depen- 
dence of signal phase fluctuation include Refs. 7, 9, 10, 14, 
and 21. It is, of course, not possible to directly measure phase 
fluctuation for noncoherent natural radio sources, as is easily 
done with coherent (e.g., spacecraft) signals. Additionally, the 
scintillation index of fluctuations from natural sources satu- 
rates of the more interesting small solar elongation angles 
as a result of the finite source size, whereas the monochro- 
matic signal phase fluctuation is not subject to such an effect 
in this region. During the period (late 1960s, early 1970s) 
when the study of natural source intensity scintillation was 
widespread, extremely valuable (monochromatic) near-Sun 
spacecraft phase fluctuation data was routinely being acquired 
by the National Aeronautics and Space Administration 
(NASA)/Jet Propulsion Laboratory (JPL) Deep Space Net- 
work for other purposes (primarily tracking system perfor- 
mance monitoring), but was not utilized in the investigation of 
solar wind density fluctuation. 

The first usages of single-frequency spacecraft data to probe 
the solar corona and solar wind consisted of a spectral broad- 
ening experiment with Mariner IV (Ref. 22), range experi- 
ments with Mariner 6 and Mariner 7 (Ref. 23), and a Mariner 6 
and Mariner 7 experiment utilizing a group-phase technique 
(Ref. 24). However, the first investigation of solar wind 
density turbulence which utilized single-frequency spacecraft 
signal phase fluctuation was an analysis of the 1975 solar 
conjunction data of the Helios 1, Pioneer 10, and Pioneer 11 
spacecraft (Refs. 25, 26). Results of this and subsequent 
investigations (Refs. 27-30) strongly suggested that signal 
phase fluctuation is radially proportional to columnar electron 
density. 

This paper presents a direct comparison of concurrent 
spacecraft signal phase fluctuation and columnar electron 
density measurements obtained during the 1976 Viking solar 
conjunction. The data argue in favor of radial proportionality 


between signal phase fluctuation and columnar electron 
density. 


II. Data 

The data presented in this paper were generated during July 
through December 1976. During this period the two Viking 
Orbiter spacecraft underwent solar conjunction with a mini- 
mum Sun-Earth-Spacecraft angle of approximately 0.25 deg 
reached on November 26, 1976. The daily rate of change of 
the Sun-Earth-Spacecraft angle was approximately 0.3 deg, 
with the Sun-Earth-Spacecraft angle at the beginning of July 
being 47 deg, and 10 deg at end of December. During this 
entire conjunction period, the NASA-JPL Deep Space Net- 
work (DSN) maintained essentially continuous coherent com- 
munications links with both spacecraft. The ground-transmitted 
(“uplink”) wavelength to the spacecraft was 13 cm (2.1 GHz, 
S-band). While the spacecraft-transmitted (“downlink”) 
wavelengths were 13 and 3.6 cm (8.4 GHz, X-band). 

During these tracking periods the DSN generates frequent 
measurements of doppler data (accumulated doppler phase) 
and less frequent measurements of spacecraft range (signal 
time delay). The doppler frequency referred to here is that 
frequency component arising from the relative motion of the 
spacecraft with respect to the ground antenna. These data are 
used for both navigation (except at small solar offset dis- 
tances) and radio science (primarily at small solar offset dis- 
tances). The doppler data are generated by comparing a very 
precise (5 X 10”*'^ stability over a day) ground reference 
frequency to the received frequency. This received frequency 
is simply the (ground) transmitted carrier frequency (S-band 
only, 2.1 GHz) which has been received, coherently multiplied 
(by 240/221 at S-band, 880/221 at X-band), and retransmitted 
by the spacecraft. 

In the course of generating these data, the DSN auto- 
matically calculated low-frequency (~10~^ Hz) S-band phase 
fluctuation information for purposes of system performance 
monitoring. In the algorithm used during this time period, a 
fixed number (usually 15) of doppler frequency samples 
(usually 60-second averages) were first differenced with the 
predicted doppler frequency which is generated by the naviga- 
tional orbit determination process. A rms doppler frequency 
fluctuation (“noise”) was then computed after fitting with 
a least squares linear curve fit to remove any possible trajec- 
tory (i.e., low-frequency) errors. 

A numerical simulation of the computational algorithm 
used (Ref. 31) provides an approximate relationship between 
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rms doppler phase fluctuation (p and rms doppler frequency 
fluctuation a^: 

0 (r) 1.67 T • (r) 

where 

T = doppler averaging time. 

Typically, a DSN ground station generates frequent doppler 
measurements over a daily period of approximately 4 to 8 
hours, referred to as a “tracking pass.” The doppler noise data 
discussed in this paper have been computed in Hz as described 
previously in this section, and then averaged over a tracking 
pass to produce a single composite value (“pass average”) for 
each daily tracking period. These Viking doppler noise data 
have already been extensively analyzed (Refs. 27, 28) and a 
more detailed explanation of the method of computation 
appears in Ref. 3 1 . 

The DSN generates precision (1.7 m) range (time delay), 
data via the modulation of a sequential binary coded signal 
onto the same carrier frequency described previously. At the 
spacecraft the range code is demodulated, and the carrier 
signal is coherently multiplied and remodulated with the range 
code and retransmitted by the spacecraft to the receiving 
station. The received S- and X-band range codes are compared 
to a reference range code which has been modulated with the 
extracted doppler frequency. The spacecraft range is extracted 
from this comparison. Appropriately multiplying and differ- 
encing the S- and X-band measurements immediately results in 
a direct measurement of columnar electron density (Ref. 32). 
Beneficially, all nonfrequency dependent effects are differ- 
enced out in this process. 

The range data are dual-frequency (S-band minus X-band) 
range delay in nanoseconds. For each tracking pass, anywhere 
from 1 to as many as 1 00 individual range measurements were 
available; a composite range delay was constructed by averag- 
ing the highest and lowest value from each tracking pass. Since 
the maximum variation in range delay measurements per track- 
ing pass was typically between 2 and 20 %, it is considered that 
the procedure was more than adequate to secure a reasonable 
daily composite range delay value. 

III. The Data Model 

Figure 1 provides the correlation between the concurrent 
range delay and doppler noise pass average data; as is seen, an 
empirically selected proportionality constant of 5 X 10“* Hz 
s“* provides a very good fit to the data. Figure 2 presents both 
range delay and doppler noise vs day of year, while Fig. 3 


presents the same data vs signal path offset distance (in solar 
radii). In both Figs. 2 and 3, the range data have been con- 
verted to Hz via the proportionality constant (5 X 10^ Hz s~*) 
shown in Fig. 1. 

The model shown in Figs. 2 and 3 represents the signal path 
integration of an electron density model N (r) =2.21 X 
108^-6 j 55 X , where the units are cm“^. The 

coefficient and the radial exponent (-2.3) of the second term 
was determined by a least squares curve fit to the phase 
Auction data in Fig. 2 (Refs. 27, 28), while the coefficient 
of the first term was determined from Viking and Helios phase 
fluctuation measurements between 2 and 6 solar radii (Fig. 4, 
from Ref. 30). The radial exponent (-6.0) for the near corona 
term derives from solar eclipse measurements (Ref. 33). The 
constant shown in Figs. 2 and 3 is equal to the speed of 
light (cm sec“*) times the columnar content (cm"^) per cm of 
electron density induced signal delay. 

The comparison of range delay to doppler noise for the 
data acquired from November 28, 1976 to December 5, 1976 
(day of year 333 to 340) appears to differ somewhat from the 
comparison over the remainder of the period and hence (this 
range data) has been marked with a different symbol in the 
figures. Possible explanations for this difference are (1) a 
region of enhanced fluctuation-to-density ratio, or less likely 
( 2 ) systematically erroneous range measurements. 

The observations here discussed are columnar quantities; 
for ease of comparison to other electron density (radial) dis- 
tribution work, one desires the generating point source elec- 
tron density distribution. To facilitate the comparison, a 
closed form approximation of the integral (of the solar wind 
component of electron density) 

( 1 ) 


where 

R = signal path 

= Earth-spacecraft distance 
r = heliocentric signal distance 

is desired. Figure 5 details the appropriate geometrical con- 
figuration. Equation (1) is rewritten as: 

7T/2+Q: 

(r^ sinay^^'''^^ I (cos co)^ rfco ( 2 ) 

''a-nl2 
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where a = Sun-Earth-Spacecraft angle, j3 = Earth-Sun- 
Spacecraft angle, = Earth-Sun distance, and: 

cu = tan“ ^ (R - cos ajr^ sin a) (3) 

A Macluarin series is utilized to expand (cos co) 

(cos co/ s 1 (-f) (-1.1 f) -t • • • (4) 

After the indicated integration one has: 

J 0 (sin F (a, 0) (5) 

where 

F(«,® a 1-0.05 

- 0.00275 j W- «/2 * - (a - »/2)’ 1 

In this work, all data fitting (Eq. 5) makes use of the two 
independent geometrical parameters a and 0, which directly 
relate to the signal path solar offset distance and signal path 
length, respectively. Much work in determining electron 
density distribution has customarily relied on usage of only 
the dominant geometrical parameter a. Significant errors can 
result in experimentally determining radial dependence via 
reliance on only one geometrical parameter (Ref 34). 

It is noteworthy that this electron density model is in good 
agreement with K-coronameter measurements at (Refs. 
35-37), pulsar time delay measurements at 10 /-g (Refs. 38, 

39) , and in situ electron density measurements at 1 AU (Ref 

40) . The model gives a significantly' higher electron density at 


lOZg than spacecraft range data from Helios 1 and 2 (Ref 41) 
and Mariners 6 and 7 (Ref 23). Table 1 presents various model 
determinations at lOr^. 

IV. Discussion 

Since at least the late 1960s, the theoretical relationship 
between columnar phase fluctuation and columnar electron 
density has been extensively hypothesized, but definitive con- 
current measurements of the two parameters have not (pre- 
viously) been available to test the various hypotheses. Recent 
efforts in theoretically analyzing the radial dependence of 
phase fluctuation (Refs. 19, 42) have produced phase fluctua- 
tion models which are proportional to the product of density 
fluctuation and solar wind (radial) velocity times a factor of 
Under the common assumptions of 
(1) proportionality between density and density fluctuation, 
and (2) the conservation of particle flux, such phase fluctua- 
tion models are immediately seen to be proportional to 
and hence are independent of the actual radial dependence of 
electron density. 

These theoretically derived phase fluctuation models of 
Woo and Callahan are not compatible with the concurrent 
measurements of Viking doppler noise (phase fluctuation) and 
range (time delay) as presented in this article. In particular, the 
models fail in the near-corona region, where the electron 
density radial dependence changes from to r~^. This can 
be seen in the near-corona region of Fig. 3 and in Fig. 4, where 
the phase fluctuation data are well fit by a model of functional 
form r"® (« / dR). Clearly, a model of the form ® 
could not reasonably fit these data. 

In way of possible explanation of this incompatibility, it is 
to be noted that analyses performed much earlier and which 
employed geometrical optics (Refs. 11, 14, and 21) predicted 
this now observed radial proportionality between signal phase 
fluctuation and columnar electron density, under the very 
reasonable assumption of a linear transverse fluctuation scale 
= Kr), as would directly be expected from a radially 
outflowing solar wind (Ref. 21). 
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Table 1. Electron density models evaluated at lOr^ 


Source 

electrons cm 

Type of 
experiment 

This article 

8000 

Dual-frequency 

range 

Ref. 41 

4500 

Helios 1 and 2 
S-band range 

Ref. 23 

5600 

Mariner 7 
S-band range 

Ref. 39 

8000 

Pulsar time 
delay 

Ref. 38 

8400 

Pulsar time 
delay 


S-BAND DOPPLER NOISE, Hz 


ORIGINAL PAGE IS 
‘ OF POOR QUALITY 



S-BAND MINUS X-BAND RANGE DELAY, NANOSECONDS 


Fig. 1. Correlation between a-^ = 13-cm signal phase fluctuation and integrated electron density. The solid data 
points are from July 2 to November 27, 1976, and December 6 to December 20, 1976, while the hollow points are from 
November 28, to D^ember 5, 1976. The solid line represents the proportionality constant 5 x 10* Hz/s, while the 
dashed line represents 1.3 x 10^ Hz/s 



S-BAND DOPPLER NOISE AND S-BAND MINUS X-BAND RANGE DELAY, Hz 
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DAY OF YEAR, 1976 

Fig. 2. Signal phase fluctuation and integrated electron density vs day of year, 1976. The solid points are = 13-cm Doppler noise in Hz, 

while hollow points are dual-frequency range in nanoseconds multiplied by the proportionality constant = 5 x 10* Hz/s. Hollow circles 
are data from July 2 to November 27, 1976, and December 6 to December 20, 1976, while hollow triangles are data from November 28 to 
Decembers, 1976. The constant K, = 3.93 x 10^'cm~^s~^ is the number of columnar electrons per second of dual-frequency range delay 
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2.21 X 10® . 1.55 X 10* 


SIGNAL CLOSEST APPROACH DISTANCE, tAq 

Fig. 3. Signal phase fluctuation and Integrated electron density vs solar offset distance. 
The solid points are = 13-cm Doppler noise In Hz, while the hollow points are dual- 
frequency range In nanoseconds multiplied by the proportionality constant = 
5x10* Hz/s. Hollow circles are data from July 2 to November 27, 1976, and December 6 to 
December 20, 1 976, while hollow triangles are data from November 28 to December 5, 1 976. 
The constant Kj = 3.93 x 10^^ cm~^s~' is the number of columnar el^rons per second of 
dual-frequency range delay 
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Fig. 4. Signal phase fluctuation In the Inner corona. 
Circles are data from the Viking 1976 solar conjunction, 
diamonds are data from the Helios 2 1976 solar conjunc- 
tion, and squares are data from the Helios 1 1975 solar 
conjunction 


Fig. 5. Sun-Earth-Spacecraft geometrical con- 
figuration. r = heliocentric distance to signal 
path, r, = Earth-Sun distance, R = Earth- 
Spacecraft distance, a = Sun-Earth-Space- 
craft angle, and = Earth-Sun- Spacecraft angle 



